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SYSTEM AND METHOD FOR CONK BEAM VOLUME COMPUTED 
TOMOGRAPHY USING CIRCLR-PLUS-MULTIPLB ARC ORBIT 
Cross-Rcfcrcncc to Related Patents and Applications 

The applicant is a named co-inventor in U.S. Patent No. 5,999,587 and the named 
■> inventor in U.S. Patent No. 6.075,836 and co-pending U.S. Paten, Application Nos. 
09/589.115 and 09/640,713. all of which concern subject matter related to the present 
invention. The disclosures of all of those patents and applications are hereby incorporated by 
reference in their entireties into the present disclosure. 
Field of the Invention 

The present invenlion is directed to a system and method for reconstruction of images 
from cone beam volume computed tomography (CBVCT) and more particularly to such a 
system and method in which the data are taken over an orbit having a circle and two or more 
aros. 

Description of Related Art 

Among all possible applications of the Radon transform, computed tomography (CT) 
applied in 2-D medical and non-destructive test imaging technology may be the one that has 
achieved the greatest success. Recognizing the demand for saving scan time in the currently 
available 2-D CT and consequently greatly improving its functionality, the implementation of 
CBVCT has been investigated for the past two decades. 

The intermediate function derived by Grangeat (P. Grangeat, "Mathematical 
Framework of Cone Beam 3D Reconstruction via the First Derivative of the Radon 
Transform," Mathematical Methods in Tomography, Lecture Notes in Mathematics 1497. G. 
T. Herman ct al, eds.. New York: Springer Verlag. 1991, pp. 66-97) establishes a bridge 
between the projection of. 3-D ohjee. and its 3-D Radon transform and is much more 
numerically tractable than pa-vim,sly known intermediate lunclions. With the prngn** in 
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understanding the so-called data sufficiency condition Tor an cxacl rcconstruclion, a few cone 
beam non-planar scanning orbits, such as dual orthogonal circles, helical, orthogonal circle- 
and-linc, non-orthogonal dual-ellipse, orthogonal circle-plus-arc, and even general vertex 
path have been proposed. Correspondingly, the analytic algorithms to exactly reconstruct a 3- 
D object based upon those non-planar scanning orbits have also been presented. 

Generally, a cone beam filtered back-projection (FBP) algorithm can make cone beam 
reconstruction much more computationally efficient and more easily implemented in a multi- 
processor parallel computing struoture. Hence, an FBP cone beam reconstruction algorithm is 
desirable in practice, and Feldkamp's algorithm (L. A. Feldkamp, L. C. Davis, and J. W. 
Kress, "Practical cone-beam algorithm," J. OpL Soc. Am. A, Vol. I, pp. 612-619, 1984) for 
the circular orbit is the earliest example. Obviously, Feldkamp's algorithm violates the data 
sufficiency condition, and an accurate reconstruction without intrinsio artifacts is available 
only in the central plane overlapping the circular orbit plane, so that some acouracy on the 
off-central planes has to be sacrificed. Although proposed independently, many algorithms of 
the prior art featured a common structure of shift variant filtering (SVF) followed by cone 
beam back-projection. Only l-D ramp filtering is employed in Feldkamp's algorithm, but a 
cascade of 2-D operations, such as weighting, 2-D projection, differentiation and 2-D back- 
projection, arc involved in the shin variant filtering. The complexity of the SVF (OfN 4 ) ) is 
higher than that of the l-D ramp filtering of Feldkamp's algorithm (OfrflogN) ). Another 
important common rcaturc possessed by many algorithms of the prior art is a normalized 
redundancy function (NRF) adopted to compensate for the multiple intersections of the 
projection plane with the source trajectory. Recently, that kind or algorithm has boon 
extended to a more general situation in which an arbitrary vertex path is involved as long as 
the data sufficiency condition is satisfied. Apparently, Ihc NRF is dala-acquisition-orbil- 
depeiulenl and has discontinuities in data acquisition orbits which moel the dala sufficiency 
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condition, bul i( can be analytically calculated for cither a specific data acquisition orbit or 
even an arbitrary vertex path. On the other hand, the algorithm by Hu (H. Hu, "A new cone 
beam reconstruction algorithm for the circlc-and-linc orbit," Proceedings of International 
Meeting on Fully 3D Image Reconstruction in Radiolog), and Nuclear Medicine, pp. 303-3 1 0, 
5 1995; and H. Hu, "Exact regional reconstruction of longitudinally-unbounded objects using 
the cirole-and-line cone beam tomographic system," Proc. SP/E, Vol. 3032, pp. 441-444, 
1997) for an orthogonal circle-plus-Iine orbit is promising in saving computation resource, 
since a window function, instead of the NRF. is employed for the cone beam reconstruction 
from the projection data acquired along the line orbit. 
"0 Due to mechanical feasibility, a circular x-ray source trajectory is still the dominant 

data acquisition geometry in all commercial 2-D/3-D CT systems cuirentiy available. Based 
upon a circular source trajectory, a number of data acquisition orbits can be implemented by 
either moving the table or tilting the CT gantry. An orthogonal circle-plus-arc orbit has been 
presented. It possesses advantages that can not be superseded by other "circle-plus" 
15 geometries, especially in the application of image guided interventional procedures requiring 
intraoperative imaging, in which the movement of a patient table is to be avoided. Further, it 
can be easily realized on a C-arm-bascd imaging system, which is being used more and more 
for tomography in recent years. The orthogonal circle-plus-an> orbit can be realized by 
acquiring one set of 2-D cone beam circle projections while rotating an x-ray source and a 
20 2D detector on a circular gantry and then acquiring another set of 2-D cone beam arc 
projections while lilting the gantry along an arc which is orthogonal to and coincident with 
the circular orbit at the same radius. The exact CBVCT reconstruction algorithm associated 
with that circle-plus-arc orbit is not in the FBP form. The rcbinning process involved in the 
algorithm requires storage for all information in the Radon space, and makes the CBVCT 
25 ^construction computationally inefficient, further, the arc sub-orbit provides information 
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covering ils Radon sub-domain only once, but Ihc circular sub-orbit provides information 
covering ils Radon sub-domain Iwioc. Thai unbalanced coverage in the Radon space may 
result in non-uniformity of noise characteristic in reconstructed images. 

A particular application of the present invention is in the detection of lung cancer and 
other malignancies. CT scanning plays a central role in much of the Ihoracic imaging used in 
detection of lung oancer and other malignancies. CT is non-invasive, easy to perform, and 
usually straightforward to interpret. It is either the primary modality or the referral modality 
for the detection of pulmonary masses (primary and metastatic), non-invasive staging of 
primary bronchogenic carcinoma, and for detection of major complications of malignancies, 
particularly pulmonary emboli, and infections. However, present helical CT has three major 
technical shortcomings. First, helical CT scans require a long or multiple breathholds for 
whole lung imaging, depending on slice thickness. Second, slice thickness vs. coverage vs. 
scan time tradeoff: programming thinner slices increases scan time or decreases coverage. 
The spatial resolution is not isotropic; through plane resolution is limited by slice thickness 
and a few times lower than that of in-plane. Third, the clinically achievable in-plane 
resolution for a large FOV, such as whole lung imaging, is limited and less than or equal to 
l.Olp/mm. 

CT of the chest is a potential screening tool for lung carcinoma. While screening 
programs based on conventional x-rays had poor sensitivity and diagnosed most carcinomas 
after the window of surgical cure had passed, CT scans reveal nodules below 1 centimeter 
with higher potential cure rates. A drawbaok of screening CT is poor specificity. Benign 
sub-centimeter nodules arc common (non-calcified granulomas, intrapulmonary lymph nodes, 
focal regions of atelectasis). The best diagnostic algorithm post-discovery of sub-centimeter 
nodules is unclear. Universal resection seems impractical. Potential diagnostic algorithms 
include evaluating the nodule enhancement, bonier characteristics, and growlh. In all of 
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these cases, accurate depiction of a small nodule is necessary. Helical CP, while readily 
delecting these nodules, has partial volume averaging problems in accurate characterization. 
11 would therefore be desirable to provide a scanning system and method with sub-millimeter 
isotropic resolution, which would potentially better characterize the density and size of these 
5 small nodules. Accurate size measurement would allow short-term follow-up to evaluate for 
growth. 

While CT screening for bronchogenic carcinoma in the high-risk population may or 
may nOt be clinically beneficial and economically practical, chest CT for the detection of 
metastases is commonly performed. CT is performed at the time of initial diagnosis, as 

10 interval monitoring for detection of disease, and as follow-up of detected nodules whioh are 
not initially resected. In all cases, improved detection and characterization, particularly that 
of interval growth, should be clinically beneficial. 

Three image intensifier (Il)-based cone beam reconstructions for volume lung 
imaging have been reported before. However, all II-based CBVCT for volume lung imaging 

15 suffers from inaccurate reconstruction due to the use of a single circle cone beam acquisition 
geometry and its corresponding approximating algorithm by Feldkamp et al, in addition to a 
limited performance of the II-CCD imaging chain. The best low contrast detectability of the 
II-based cone beam CT for volume lung imaging is 10 HUs for a 3 mm object. 

20 
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Summitry of llic Invention 

It will be readily apparent from the above that a need exists in the art to overcome the 
above-rioted limitations of the prior art. II is therefore an object of the invention to satisfy the 
data sufficiency condition while achieving a more balanced coverage. It is another object of 
5 the invention to do so in a computationally efficient manner which can be adapted to parallel 
cone beam reconstruction. 

To achieve the above and other objects, the present invention is directed to a system 
and method for reconstructing images from data taken over a circle and two or more arcs. An 
FBP reconstruction algorithm is presented for reconstructing the images. 

10 The efficiency of reconstruction is critical for the application of CBVCT in the image- 

guided interventional procedures, and the reconstructed images with uniform noise 
characteristic are desired in practice. In order to overcome the previously mentioned 
shortcomings of the circle-plus-arc orbit and its associated Radon Transform-based 
reconstruction algorithm, a cirole-plus-two-arc orbit and an analytic FBP cone bean 

15 reconstruction algorithm are used. The result given by Hu for the circular cone beam 
projections is directly incorporated. For the cone beam projections acquired along the arc 
orbits (namely, arc cone beam projections), originating from the equation established by 
Grangeat and the inverse Radon transform, an analytic reconstruction solution is obtained. 
That solution is different from known solutions because a window function, instead of an 

20 NRF, is employed to compensate for the multiple intersections of the projeotion plane with 
the x-ray source trajectory. Since its support in the Radon domain is very limited, the window 
function of the present invention significantly reduces the computational cost of the 
reconstruction from the arc CD projections. 

Most objects to be reconstructed in medical or non-deslruclivc x-ray CT are 

25 longitudinally unbounded. I lence, a cone beam reconstruction algorithm should address such 



6 



WO 02/30282 „„ „ 

PCT/DS01/32011 

a truncation problem. In order lo solve (he so-called truncated cone beam projection, several 
methods have been proposed. It has been demonstrated that a finite region of interest (ROI), 
for which the extended data sufficiency condition is satisfied, can be reconstructed 
accurately, although that finite ROI is slightly smaller than the ROI which can be scanned by 
5 a detector. The circlc-p| us -lwo-arc orbit and its associated cone beam FBI> reconstruction 
algorithm in the present invention are intrinsically capable or dealing with the truncation 
problem, and its thorough evaluation is accomplished herein. 

The circle-plus-arcs orbit possesses advantages over other "circle-plus" orbits for the 
application of x-ray CBVCT in image-guided interventional procedures requiring 
10 intraoperative imaging, in which movement of the patient table is to be avoided. A cone beam 
circle-plus-two-arc orbit satisfying the data sufficiency condition and a filtered back- 
projection (FBP) algorithm to reconstruct longitudinally unbounded objects is presented here. 
In the circle sub-orbit, the algorithm employs Feldkamp's formula and another FBP 
implementation. In the arc sub-orbits, an FBP solution is obtained originating from 
15 Grangeat's formula, and the reconstruction compulation is significantly reduced using a 
window function to exclude redundancy in Radon domain. The algorithm's merits include the 
following: Only 1-D filtering is implemented even in a 3-D reconstruction, only separable 2- 
D interpolation is required to accomplish the 3-D back projection, and the algorithm structure 
is appropriate for parallel computation. 

The present invention has the following characteristics and advantages. A fiat panel 
detector (FPD) can be used. The invention can incorporate scattering correction and volume 
of-inlerest (VOI) reconstruction. The present invention can be used for medical imaging, 
nondestructive testing or any other purpose in which such imaging is desired. 

In the reconstruction algorithm of the preferred embodiment, all the components are 
25 in a filtered backprojedion format. That reconstruction algorithm is more computationally * 
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efficient Hum those of the prior art and is ready for parallel cone beam reconstruction. Thai 
algorithm can be used lo provide an exact reconstruction of a longitudinally unbounded 
object. The CBVCT reconstruction of the preferred embodiment is the 3D matrix or 
attenuation coefficient distribution of a 3D object. 

In the present invention, the data arc taken through a scan such as a quasi-spiral scan. 
To achieve the fastest scan, a simplified scan, such as only tilt in plus circle scan, can be used 
to satisfy the data sufficiency condition. The second set of arc projection scans (gantry tilt-out 
scans) is optional to improve image quality. The total acquisition time can be reduced by 
decreasing the sampling rate on the arcs or by using only a gantry-tilt-in plus a circle scan 
during the quasi-spiral scan. 

The present invention offers the following particular advantages when used to detect 
lung cancer. First, the present invention requires a much shorter volume scanning time 
relative to helical CT. In a single volume scan, an entire acquisition can be performed. The 
present invention can improve acquisition efficiency by a factor of 25 for 1 mm slice 
thickness per volume scan vs. a single ring helical CT. Assuming a 25 cm segment to be 
scanned for a whole lung imaging and I mm/slice, the present invention can be at least 24 
times faster than a single ring detector helical CT and 3 (for gantries with 0.5 sec/revolution) 
to 6 times faster than a multi-ring detector helical CT. The fast volume scan eliminates the 
respiratory misregistration problems, such as those caused by the requirement that the patient 
hold his or her breath, and is less sensitive to patient motion. 

Second, the present invention can provide isotropic resolution in the x, y and z 
directions and provide true 3D reconstruction images. The spatial resolution of PI'D-bosod 
CBVCT is limited by the fineness of our detector array, not by collimation. An PPD-boscd 
CBVCT achieves spatial resolution on the order of 1-2 lp/mm in routine mode. The present 
invention can provide higher resolution in all three directions than a helical CT. 
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Third, the embodiment with ultra-high resolution VOI reconstruction can provide true 
3D tomographic reconstruction with spatial resolution approaching that of screen-film 
projection imaging, but with 50 -100 times better contrast resolution than projection imaging. 
This spatial resolution capability cannot be achieved in any current helical CT. 
5 In addition, the present invention can more efficiently use x-ray tube output and 

greatly reduce the tube loading requirement. This will reduce the manufacture cost of CT 
tubes because a very heavy duty and very costly x-ray CT lube ($60,000 - $100,000/tube) 
may not be needed, and/or the operating cost because the life of a CT tube will be many times 
longer. 

10 The present invention thus improves the sensitivity and specificity of lung cancer 

deteotion as well as other types of cancer detection. In addition, it will highly significant to 
the early detection and management, not only of lung cancer, but also of other malignancies. 

There are several radiological or biologioal characteristics of carcinoma that can be 
imaged. First, carcinoma has different x-ray linear attenuation coefficients from surrounding 

15 tissues. Second, carcinoma has a substantially higher volume growth rate compared to a 
benign tumor, which lacks growth. Third, carcinoma has border patterns distinguishable 
from those of a benign tumor. Fourth, benign tumors show different contrast enhancement 
after intravenous contrast injection. Fifth, the presence of neovascularity can indicate cancer. 
Conventional cancer detection techniques such as chest projection imaging rely mainly on the 
20 first characteristic and partially use the third characteristic for cancer detection. Since 
mammography is a two-dimensional static imaging technique, it cannot provide any 
information regarding characteristics 2, 4, or 5. The present invention, by allowing fast scans 
and permitting the use of contrast injection if desired, can be used to detect cancers in 
accordance with all five characteristics. 
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CT scanning is n key modality for delecting pulmonary malignancies. It can detect 
lesions as small 2-mm diameter. It is, however, imperfect for detection of nodules for the 
following reasons: 

Nodules may not be imaged if the lungs cannot be scanned in a single breathhold. 

5 Respiratory misregistration occurs when a CT scan of the lungs is acquired in several 
different breathholds. Because patients do not reliably hold their breath in the same phase of 
respiration, and because pulmonary lesions move cranially or caudally with respiration, a CT 
scan composed of slices obtained from different breathholds may fail to detect a lesion 
because that lesion was never imaged. The present invention permits scanning the entire 

10 lungs in a single breathhold and thus can eliminate this source of detection error. 

Nodules may be present on the CT images but fail to be recognized by the interpreting 
radiologist. A retrospective review of nine patients with missed lung cancer on CT found five 
missed tumors that were peripheral and <3mm in diameter and four central tumors measuring 
up to 8 mm in diameter. These small peripheral nodules were likely not seen while the larger 

IS central nodules were not recognized set against the background of the larger complex 
branching vessels. Review of a CT dataset electronically, and in planes other than the axial 
plane might also prove to have further increase in sensitivity for nodule detection. The 
present invention provides the first system capable of scanning the entire chest with sub- 
millimeter isotropic resolution. Isotropic resolution with sub-millimeter resolution in all 

20 directions would be ideally suited for electronic interpretation in axial, oblique, coronal and 
sagittal planes. 

Partial volume averaging with adjacent lung can make small pulmonary nodules 
difficult or impossible to detect by helical CT. Helical CT of canine metastatic ostcosaroma 
found 44% of metastases <5 mm vs. 91% of metastases >5mm. Usually, helical CT 
2? reconstructs images at an interval approximately equal to the eollimation, 5-7 mm. Some 



10 



W0 02/30282 PCT/US01/32011 
clinically relevant nodules ore smaller than (he sliee thickness. Reconstructing images at a 
smaller reconstruction interval increases the sensitivity for lung nodule detection. This is due 
to the non-linear slice sensitivity profile of helical CT reconstruction. These overlapping 
reconstructions have a belter chance of placing small nodules in the center of the slice where 
they will be displayed with higher density and be more easily seen. The present invention can 
overcome this problem because slices at <1 mm thick would essentially eliminate partial 
volume averaging and also assure that a nodule larger than 3mm would have a slice through 
its approximate center. 

Nodule size is difficult to measure accurately by helical CT. The apparent size of a 
pulmonary nodule depends on the thickness of the slice and where the slice is reconstructed 
relative to the nodule. Accurate size measurements of the nodules are necessary to detect 
small amounts of growth in short-term follow-ups. A 3 mm diameter nodule growing to 4 
mm diameter has more than doubled in volume. Because the detection of small nodules is 
becoming increasingly common due to helical CT, it is likely that imaging algorithms will 
need to incorporate follow-up of small nodules for growth. Accurate sizing would be 
essential. CBVCT will provide 0.125 - 0.7 mm voxel size and will allow accurate 
measurements of nodule size and nodule volume. 

Small nodule density (attenuation coefficient) is difficult to measure accurately by 
helical CT. The apparent density of a pulmonary nodule in helical CT depends on the 
position of the nodule relative to the position of the reconstructed slice. The relative 
movement of the slice by one or two mm may make a calcified nodule appear non-calcified. 
For nodules smaller than the slice thickness (routinely 5-10 mm in helical CT), there is partial 
volume averaging of the nodule with adjacent air and an accurate density can not be 
determined. The nodule density is useful for characterization in two major respects. One is 
the detection of calcification indicating benignity. The second is that malignant pulmonary 
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nodules appear to have more rapid contrast enhancement than benign nodules. Sub- 
millimctcr thick slices, achieved by CBVCT, will allow accurate density measurements of 
small nodules without partial volume averaging and without necessity for post-processed 
overlapping reconstructions. This should better detect calcification, and more accurately 
5 characterize the amount of enhancement. 

Fine spiculalions and other nodule border characteristics arc best determined with 
high resolution CT. On helical CT scanners, this requires locating the nodule prospectively, 
as it is impractical to acquire high-resolution CT I -mm thick slices throughout the entire 
lungs. CBVCT would acquire high-resolution images through every nodule without prior 

10 knowledge of its location or the need for technologists or physician localization during the 
scan. The ultra-high resolution VOI reconstruction mode of CBVCT will provide even 
higher resolution for target imaging after the survey lung imaging of CBVCT with lower 
resolution. This mode may be even more useful for characterizing nodule border. The value 
of universal high resolution CT for characterizing benign vs. malignant nodules may also 

15 prove beneficial. 

A particular implementation of CBVCT provides high contrast resolution of 0.7 - 4 
lp/mm, and low contrast detectability of 3-5 CT number within a short breath hold (2-8 
seconds). Such an implementation preferably includes an appropriate 2D deteotor system 
which has a high detection quantum efficiency (DQE), high dynamic range, high spatial 

20 resolution, minimal geometric distortion, and which is capable of high image acquisition rates 
with little image lag and excellent linearity, ll also preferably includes a data acquisition 
scheme that will result in a complete set of projection data with little additional mechanical 
complexity. This provides an exact cone beam reconstruction algorithm which is based on 
the Complete set of data, thereby permitting imaging in a large l-OV (for example 14" - 16"). 
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A third aspect which is preferably included is x-fuy scatter control and correction techniques 
to further improve low contrast dctcclabilily. 

Tlic present invention expands the application of CBVCT from angiography to 
volume lung imaging and other applications that require soft tissue differentiation. CBVCT 
can potentially be applied to pulmonary emboli detection, liver cancer detection, volumetric 
brain perfusion, diagnosis of acute stroke, and colon cancer detection, etc. 

For lung cancer and other malignancies, the present invention has application to 
malignancy detection, monitoring, management and treatment and in particular to the 
development of treatment plans. 
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Brief Description ofthc Drawings 

A preferred embodiment of the present invention wilt be set forth in detail with 
reference to the drawings, in which: 

Pig. I is a schematic drawing showing a cone-beam projection; 

Fig. 2 is a schematic drawing showing a Radon plane; 

Fig. 3 is a schematic drawing in Radon space showing the inability of the circular 
orbit alone to satisfy the data sufficiency condition; 

Fig. 4A is a schematic diagram showing the ability of the circle-plus-two-arc orbit to 
satisfy the data sufficiency condition; 

Fig. 4B is a diagram of a quasi-spiral scan used to implement the circle-plus-two-arc 

orbit; 

Fig. 4C is a diagram of a scan used to implement a circle-plus-two-line orbit; 

Fig. 4D is a diagram of a scan used to implement a helical orbit; 

Fig. 5 is a schematic diagram showing the coordinate system and parameters used in 
the reconstruction of the circular projection; 

Fig. 6 is a schematic diagram showing the coordinate system and parameters used in 
the reconstruction of the arc projections; 

Fig. 7 is a plan diagram of a system on which the preferred embodiment is 
implemented; 

Figs. 8-1 1 show stages in the operation of the system of Fig. 7; and 

Figs. 12A and 12B show a setup for taking scout images for scatter correction. 
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Detailed Description of (lie Preferred Embodiment 

A prcrerred embodiment of the present invention will now be set forth in detail with 
reference to the drawings. First, the reconstruction algorithm will be derived. Second, a 
device on which the reconstruction algorithm can be implemented will be shown. 

The cone beam projection of a 3-D object is schematically illustrated in Fig. I, where 
O is the origin of the coordinate of the real 3-D space R\ upon which the algorithm is 
derived. Y and 2 are the local axes in the plane of the virtual detector. OS = O is the vector 
which represents the cone beam focal point S, and point A ' is the projection on the detector 
plane of A, which is a point within the 3-D object to be reconstructed, along the unit 
directional vector 



\SA\ 



The vector from O to A is r . The cone beam projection of the 3-D object f{r) is defined as: 

gfa* )«£/(*+ '<?)* (2) 
The reconstruction algorithm is derived for a longitudinally bounded object first, and its 
capability of regionally reconstructing a longitudinally unbounded object will be analyzed in 
detail later. Since most objects to be reconstructed in medical and non-destructive test 
tomography arc cylinder-like, the 3-D object to be reconstructed is assumed to be a cylinder 
whose half height is represented by h s and radius by R, 

The definition of a 3-D Radon plane P(/?,p) is schematically illustrated in Fig. 2, 

where 

#ie(siiiOcas<p t sin08inq) v cosO) (3j 
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is Ihc normal vector, and p the distance away from the coordinate origin O. In 3-D 
Cartesian coordinates, any plane can be uniquely identified by a normal vector and a distance 
away from O and is the set of all points for which 



Recognizing the existence of several versions of the data sufficiency condition, (lie 
preferred embodiment uses the one which can be the most simply expressed: All planes 
passing through the object to be reconstructed must contain a point in the scanning orbit. 
Obviously, a circular orbit, which is the simplest in practice for 3-D reconstruction, violates 
the data sufficiency condition in a way demonstrated in Fig. 3, i.e., the Radon transform of 
the circular CB projections provide no coverage on the shadowed sub-domain in the Radon 
space. In the perspective of inverse Radon transform, the sub-domain missed in the Radon 
space by the circular orbit (namely missed Radon sub-domain) has to be covered by 
additional non-circular orbits. 

The requirement for the circle-plus-one-arc orbit to meet the data sufficiency 
condition and cover the 3-D cylinder object to be reconstructed completely is known in the 
art to be (see Fig. 4A): 



where R is the radius of the cylinder object to be reconstructed, D is the radius of the arc 
orbits, and X lt1ih A is the minimum arc spanning angle range of the whole single arc. 

Since an uneven sampling is achieved by the circle-plus-arc orbit in the missed Radon 
sub-domain, some artifacts may occur in reconstructed images. Further, the arc sub-orbit is in 
odd-symmetry in the plane determined by itself. Such an odd-symmetry makes the non- 
unilbrmily of the sampling in Radon space worse, and may result in more reconstruction 



r • ft - p = 0. 



(5) 




(6) 
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artifacts. On the other hand, the redundancy function equals I within the missed sub-domain 
in the Ration space, while the redundancy Amotion equals 2 within the sub-domain covered 
by the circular data acquisition orbit. In the perspective of signal processing, that kind of 
difference in the redundancy function results in an uneven noise characteristic in 
reconstructed images. To avoid the artifacts caused by the odd-symmetry and maintain 
identical Radon information redundancy between the sub-domains covered by the circular 
and arc orbits respectively, a circle-plus-two-arc orbit is used in the preferred embodiment, 
although the invention could be adapted for a circle plus more than two arc orbits. As 
schematically shown in Fig. 4A, the circle-plus-two-arc orbit consists of a circle and a pair of 
arcs. One is called arc sub-brbit I and is represented by the solid curves; the other is called 
arc sub-orbit 2 and is represented by the dashed curves. The arc sub-orbit plane is 
perpendicular to the circular sub-orbit plane, and they are concentric at point O with the same 
radius D. It is noticed that an even-symmetry is achieved by integrating arc sub-orbit 1 and 
arc sub-orbil 2 in the arc sub-orbit plane (x 0> y 0 ) . 

In the fixed coordinate system (x, 1 y 01 zj illustrated in Fig. 4A, the circle-plus-two- 
arc orbit can be analytically depicted as 

fo c (\) = (D cos X, 0, -DsinX) 

1 XeA c =[0, 2k] (7) 
<t> 0 ,(M = (£cosX t DsinX, 0) 

| <L,(*) = (0cosX, DsinX, 0) 

l*eA tf] =|0, X mllll/ |Ul7c, x+X^J (9) 
where ^.(X) represents the circle sub-orbit, and f lf (X) the are sub-orbit I, and (|) tl (X) the 
arc sub-orbit 2, respectively. Theoretically, the minimum cone angle for each arc in the 
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eifcle-pliis-two-nrc orbit is cut clown lo half that required hy (he single ai*c in the cirelc-plus- 



one-arc orbit 



(10) 



Accordingly, the reconstruction algorithm can be most broadly expressed as 

ftr)~fr(rhf.{?) (ii) 

where f ( (r) is the component reconstructed from the sub-domain in the Radon domain 
support corresponding to the circular cone beam projections, and /„(r) the component of the 
arc cone beam projections. 

The coordinate system on which the reconstruction algorithm for the circular cone 
beam projections is derived is illustrated in Fig. 5. (x K ,y t i z t ) is the local coordinate system 
which rotates rigidly in phase with the detector. F represents a vector that determines a point 
A within the 3-D object to be reconstructed, and (y,Z) is the coordinate of the projection of 
A in the detector coordinate system. The circular sub-orbit is within the plane determined by 
( z o> x o)> and f 0 {Y,Z) >s the cone beam projection with the source focal point at O, 

It has been shown that f c (?) can be further split into: 

/c(?)-/c,(r)+/ Cl (r) (12) 

/«, (?) corresponds to the FDK algorithm and can be obtained using the following formulae 
that are modified to match the coordinate system shown in Fig. 5: 

4n W-r-x L ) 



IJ-I -.\, /)..,- ..v ( 



(14) 



/C(r.z)- j*hi u (z-z)i\(r.z) ( | 5) 



IK 



*"»™™ NHMW, 

j|a)|c/cDCxp(y'coZ) (17) 

where co 0 is the integral limit, which is determined by the spatial Sampling frequency of the 
dcleotor, in the Fourier domain. 

On the other hand, / Cj (r) can be obtained using the following formulae that are also 
modified to match the coordinate system shown in Fig. 5: 

K(r)-jp".(y)- Kfr-^.CM' (20) 
<u00- (2d 

Old 

h j*(y)~ Jycoexp(/'o)J'Vo) ^23) 

where L. is the integral limit along (he Z direction, and <o 0 is the same as that in (17). 

Both (13) and (18) are in the FBP form, and reconstruction from the circular cone 
beam projections is computationally efficient because only l-D filters //,,„(*) and /;„„(,') are 
involved in the filtering process. In Iho digital signal processing point of view. 
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''»>)=• 



d* .„ " = 0 

with (24) 



and 



i-(-i)" i ' «*o 

o„ J ..J 



c-u-i (25) 



/7 

In order to reduce numerical artifact and constrain noise, a Hamming window is implemented 
while filtering. 

It is known that the Radon transform of the circular CB projections fulfills only a 
torus in the 3-D Radon domain. As argued by Hu, the assumption that the redundancy 
function equals 2 is valid only for the Radon domain point inside the torus. With respect to 
the Radon domain point on the boundary of the torus, which corresponds to the tangential 
intersection of the projection plane with the circular orbit in the spatial domain, the 
redundancy function equals 1. Hence, the reconstruction implemented using Feldkamp's 
algorithm takes only the contribution from the Radon domain point inside the torus. In order 
to take the contribution from the Radon domain point on the boundary of the torus into 
account, the complementary term / rj (r) should be incorporated into the algorithm to 
implement the reconstruction from circular cone beam projections. 

As eluoidated above, employing the arc sub-orbits provides information in the Radon 
domain to cover the missed Radon sub-domain. Considering computational efficiency, a 
reconstruction algorithm for the arc cone beam projections in the FBP form is desired in 
practice. 

Before the transform itself is presented, some variables will be defined with reference 
to Mg. 0. The Radon plane containing S and A intersects the detecior to define a line segmenl 
having end points l)\ and D± The point of closest approach of that line segment to O is at a 
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point C The line segment connecting O to C has u length / and defines an angle <$ with the 
)>/. axis. The origin S of the cone beam is along one of the arc orbits at an angle p from the x„ 



axis. 



as: 



The 3-D Radon transform and its inverse of the object /(/") are defined respectively 



*(".P)- fyW'^-pK (26) 

f(?) =-^t l^^ sinG |r^P)l P ^ " (27) 

Based upon the geometry shown in Fig. 6 and originating from the equation established by 
Grangeat and the inverse Radon transform in Equation (27), the reconstruction algorithm for 
the arc CB projections can be written as (sec APPENDIX A) 

/«(^) = 0.5-/ O) (r)4.0.5./ aj (r) {28) 
where the factor 0.5 is to compensate for the data redundancy resulting from the double 
coverage on the missed sub-domain in the Radon domain by the two arc sub-orbits, and 
/„, (F) ( ' " (h 2} ) can be expressed in the FBP form: 



4A 



(29) 



P^p,/ ) e)=^ M <P,/,o)cos0[f| + ^l^J /((P ,,,e ) (30) 

K, ). 

/ ( CfU,0)= J J/' p .,O'^)5()'sin0 + ^cose-/>/}V« (si) 



>^) = r.,. (7r ^— and Zfrj-r-r, 



(32) 



where /~ represents a vector that determines a pninl A within the 3-1) object to he 
reconstructed, ()\/) is- the coordinate of the projection of A in the detector coordinate 
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system, )/and Z, are (lie integral limits along J'and /faxes respectively, P^()\X) ( /»//, 2/ ) 
is the cone beam projection at angle along the arc orbits, and 



1 /sin p > D (l-cos0cosP) 

I /sinf* <-D.(]+-cos0cosfl) (33) 
0 elsewhere 



is the window function derived in Appendix B to resolve the data redundancy and constrain 
5 the back-projection for the arc cone beam projections. The support of the window function 
in the sinogram domain is very limited, and the computational resources for the 
reconstruction from the arc cone beam projections can be saved substantially. Both the l sl 
and 2 nd derivative of the sinogram along the distance direction are obtained using the 1-D 
linear digital operator 

10 Notice that the algorithm structure of Equations (29>(33) look similar to the 

algorithm presented by Hu. However, there are important differences between the derivation 
of Equations (29)-(33) and that in the prior art. First, each source point defines a sphere in the 
Radon domain (namely the Radon sphere), and the diameter of a Radon sphere is determined 
by the distance between the source point and the origin of the coordinate system. The 

15 diameters of the Radon spheres along the arc sub-orbits in that algorithm are constant, but 
those along the line sub-orbit are variable. With respect to the FBP CBVCT reconstruction, a 
scries of Radon spheres with identical diameter will sample the missed Radon sub-domain 
more uniformly than a scries of Radon spheres with varying diameters. It is possible for a 
more uniform sampling in the missed Radon sub-domain to create less artifacts in the FBP 
20 cone beam reconstruction. Second, the window function (33) is distinct from that in the prior 
ad. 

The capability of regionally reconstructing a longitudinally unbounded object is 
essential for the application olTHV(T in medical or non-destructive (est imuging, since most 
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objects to be reconstructed in practice are longitudinally unbounded (that is also called the 
truncation problem). The circlc-plus-two-arc orbit satisfies the extended dala sufficiency 
condition proposed by Kudo and Saito (H. Kudo and T. Saito, "An extended completeness 
condition for exact cone-beam reconstruction and its application" IEEE Con/ Rec. 1 994 
Nuclear Science and Medical Imaging Symposium, Norfolk, Virginia, pp. 1710-1714, 1995). 
Consequently, its associated cone beam FBP reconstruction algorithm presented above 
addresses the truncation problem by employing the window function m{(3,/,0), even though 
the object to be reconstructed is assumed longitudinally bounded in its derivation. That 
means that an ROI within a longitudinally unbounded object can be exactly reconstructed if 
the ROI is smaller than the region that can be completely covered by the x-ray tube-detector 
during a scan along the circle-plus-two-arc orbit. 

On the other hand, both h M (n) and h^{n) involved in die reconstruction algorithm 
further lessen the ROI that can be exactly reconstructed. As shown in (15), the filtering by 
/^(n) is implemented latitudinally in obtaining /Jir). Since the object to be reconstructed 
is latitudinally bounded, ^,(77) incurs no contamination to /J/"). However, the filtering by 
h Ja (n) is implemented longitudinally in reconstructing / ej (r) (20), and incurs contamination 
t0 A 00 because of the longitudinal truncation. Similarly, the filtering by h Ja {n) incurs 
contamination to fjr) ( / = {I, 2} ). Fortunately, both A M (*) and A» arc square 
summablc and drop dramatically, and make the contamination depth resulting from them to 
y;,0")> A, 0 s ) and f, u (r) very limited. 

Theoretically, data redundancy can be used to improve the signal to noise ratio (SNR) 
ofn reconstructed image in CBVCT. However, unlike nuclear medicine, acquiring redundant 
projection data in an x-ray CMVCT may result in unneees.st.ry radiation to a patient. 
Therefore, a candidate scanning orbit for application in an x-ray C'HVCT should keep the 
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dulu redundancy us low as reasonably achievable (ALAKA criterion) while sutislyiiii* the data 
sufficiency condition and maintaining the image quality of a reconstructed image clinically 
acceptable. The circlc-plus-two-arc orbit with the cone beam PUP reconstruction algorithm 
presented here is one that meets the ALARA criterion. Hence, the evaluation of its 
5 performance, such as the quality of the reconstructed image as a function of arc orbit angle 
sampling interval, arc orbit spanning range, and x-ray source quantum noise levels, as well as 
its capability to regionally reconstruct a longitudinally unbounded object, is practically 
important. In order to avoid the transition between the circular data acquisition and the arc 
data acquisition, the circle-plus-two-arc orbit can be implemented through a "quasi" spiral 

10 scan. In that scan, the x-ray tube- detector mounted on a circular gantry continuously rotate. 
The circular sub-orbit is realized by acquiring 2-D cone beam projections at evenly 
distributed angular positions along one circle of the x-ray source trajectory while the tilting 
angle of the gantry is 0* (see Fig. 4B). The arc sub-orbits are realized by acquiring 2-D cone 
beam projections at both the top and bottom arcs . The quality of the reconstructed images is 

15 still acceptable while the arc sub-orbit sampling interval is only half the circular sub-orbit 
sampling interval. That means that the total turns of the "quasi" spiral scan in the circle-plus- 
two-arc orbit can be decreased significantly. Hence, the data acquisition time along the arc 
sub-orbits can be reduced remarkably. Such a significant decrease in data acquisition time is 
practically important in the application of CBVCT in the image-guided interventional 

20 procedures. 

Th capability of the cone beam FBP algorithm to regionally reconstruct a 
longitudinally unbounded object has been verified. 'Hie survival of the algorithm from the 
truncation problem is essential for its application in CBVCT. On the other hand, in the case 
of shortened arc sub-orbits that violate the data sufficiency condition, a regional exael 
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reconstruction can still be obtuined. Thai menus that the spanning range of (he are sub-orbit 
can be lessened if only an ROI within the objeel is to be reconstructed. 

In implementing the algorithm on a computer system, the reconstruction load is 
divided into several parts and run in parallel on a RACIi parallel computation system which 
is a scalable mulli-proccssor-bascd system and is provided by Mercury Computer Systems. 
Initially, a RACE with 8 upgraded processors will be used, so that the reconstruction time of 
the algorithm will be 10 - 12 minutes. Further reducing the reconstruction time by parallel 
computation to 2 minutes for 512 3 matrix reconstructions, for 288 projections with 512 x 512 
pixels per projection, can be achieved using a RACE having 16-32 processors with 1024 
Mbytes RAM at a relatively low cost. 

In a standard CT t a 3-D reconstruction is obtained by stacking a series of slices. In an 
CBVCT, a direct reconstruction of an object can be obtained. Referring now to FIG. 7, it is 
shown how a CBVCT system 700 of the present invention can be used to obtain a direct 3-D 
reconstruction of an object. It should be understood that the CBVCT scanning apparatus 700 
is illustrated in a simplified block diagram form. The invention may preferably be employed 
in conjunotion with such a CBVCT scanning apparatus to generate a 3-D reconstruction 
matrix of the object. Based on the 3-D reconstruction matrix, the desired three-dimensional 
display can be obtained. 

A C13V CT scanning apparatus examines a body P using a cone shaped radiation 
beam 704 which traverses a set of paths across the body. As shown in FIG. 7, an x-ray source 
710 and a 2-D detector 711 such as a flat-panel detector are mounted on a gantry frame 702 
which rotates around the body I* being examined. The operating voltage for the x-ray source 
is obtained from a conventional high-voltage generator 708 in such a manner that the x-ray 
source 710 produces the desired cone-shaped beam of radiation when the high-voltage is 
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applied 10 it. The high-voltage generator 708 is energized by mcuns of n power source 718. 
through a switch 716. 

A first motor 712 is also powered by the power source 718 such that it drives the 
gantry frame 702 in its orbit about the body, for example, in a clockwise direction as shown 
5 by the arrows adjacent to the frame. The power source 718 is turned on by means of switch 
720 or other conventional control devices, in order to initiate a measurement sequence. A 
speed control circuit 714 is used to control the speed of rotation of the gantry frame 702 and 
to provide an output control signal which indicates when the speed of the motor 7 12 is at the 
desired level for taking measurements. The output from the rotational control 714 may also 

10 be utilized to operate the switch 716 such that the high-voltage generator 708 is only turned 
on when the gantry frame 702 is driven at the desired speed for making measurements. 

In order to obtain the arc measurements as previously discussed, a tilt control 715 is 
utilized to cause the gantry frame 702 to tilt by a relatively small angle of ±15° to ±30°, by 
means of the gantry frame tilt motor 713. That tilting allows the acquisition of arc projection 

IS data on the perpendicular arc. Such geometry results in a complete set of data for an object 
with a 25-40 cm diameter corresponding to a 37-60 cm field at the detector 71 1 with a 
magnification of 1.5. Although the tilting of the gantry 702 is generally available in a 
standard CT gantry, to acquire arc projections, the minimal modification of a standard CT 
gantry has to be made such that the tilting of the gantry, the x-ray exposure timing and the 

20 projection acquisition are synchronized by a system control computer 724 having a clock 
722- 

Hie gantry can be based on modifications of existing equipment made by Such 
companies as Oli, Siemens, Toshiba and Marconi. Such modifications include replacing the 
one-dimensional detector with an II-CCI) detector or a silicon or selenium thin film transistor 
25 array I- PI > and the old computer system and its control interlace hoards with a new host 

2f> 
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computer and now interitioe hoards. As explained in the co-pending npplicnlions cilcd above, 
a slip ring is preferably used to permit communication between equipment on tbc gantry and 
equipment off the gantry. Initially, volume scanning speed will be only limited by the 
maximum frame rate or the real time FPD. Currently available real time FPDs have a frame 
5 rate of 15-120 frames/sec. The flat panel researchers predict that the future frame rate can be 
up to 120 framcs/sec. (IK x IK pixels/frame) and 480 frames/sec with reduced vertical 
readout lines (256 x 1 K pixels/frame). When the frame rate of the detector is increased to 480 
frames/sec. for a large size FPD in the future, the volume scanning time of entire lungs will 
shorten to 1-2 seconds depending on the required resolution and/or the projection number can 
10 be increased to improve image quality. Compared to II-based VTDA systems, the FPD-based 
CBVCT system represents a significant technologic advancement due to using flat plane 
detector, slip ring technology, and cone beam reconstruction algorithms Uiat result in accurate 
reconstruction. In addition, the CBVCT system can incorporate a scaleabie multiprocessor- 
based parallel computing system (8-32 processors) provided by Mercury Computer systems. 
15 A specific scanning protocol will now be described whioh implements -15° to +15° 

tilting to obtain 25 cm coverage in the z direction. This protocol consists of four steps: 1) 
Positioning gantry - Before starting the scan, the gantry tills to -15° to prepare for CPA scan; 
2) Arc-Projection Acquisition (Gantry lilt-in) - While the gantry is tilting from -15° to 0°, the 
x-ray tube and delcotor rotate, taking projections only at 0° (on the upper arc) and 18Q° (on 
20 the lower arc) of the rotation angle positions to obtain two arc projections per rotation; 3) 
Circle Projection Acquisition - When the gantry tills to a 0° tilling angle, the gantry stops 
tilting, and the x-ray tube and detector rotate to acquire multiple circle projections; and 4) 
Arc-Projeelion Acquisition (Optional Gantry lilt-out) ~ ir necessary, after completing circle 
scan, the gantry lilts from 0° to H 15°, while (he x-ray tube-detecK.r rotates, taking arc 
25 projections as in step 3. Figure 4H shows eirele-plus-are scan with six arc projections taken 
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ul the positions labeled I through 6 along the upper and lower ares. Fig. 10 shows exposure 
with the gantry tilled, while Fig. 1 1 shows exposure with the gantry not tilted. 

To reduce circle-plus-arcs CBVCT scan time , the quasi-spiral scan mode or the 
gantry is used because during the scan, the x-ray tube and detector continue rotating on the 
5 gantry while the gantry is tilling and the gantry stops lilting at 0° tilting angle lo acquire circle 
projections. The quasi-spiral scan mode eliminates the need to stop the rotation of the x-ray 
tube and the detector during the scan and reduces the transition lime between arc acquisition 
and circle acquisition. In addition, a complete set of cone beam projection data can be 
achieved using two opposite half arcs (1-4 arc projections) and a single circle scan orbit. For 

10 example, as shown in Fig. 4B, arc projections can be taken only at locations 1-4, or only at 
locations 5 and 6, corresponding to gantry tilt-in without tilt-out. With two complele arcs, 
e.g., projections at all of locations 1-6 of Fig. 4B, image quality is better. Therefore, gantry { 
tilt-out is an optional mode which can be eliminated in the interest of time constraints. In 
other words, if the imaging task requires high temporal resolution to reduce motion artifaots 

15 or to obtain dynamic information, only a gantry-tilt-in arc scan and a circle scan are required, 
which reduce the arc scanning time to half. 

In detection of lung cancer, since only 25 cm of the trunk of the body will be viewed 
per scan in the z direction, the gantry needs to be lilted ±15° at most. The volume scan time 
should bo 4-8 seconds, depending on the achievable tilt speed, how large the segment in the z 

20 direction is actually viewed and the acquisition rate of the dotcctor. The system provides 
computer-controlled gantry tilt and synchronized x-ray exposures with 2 cxposures/scc for 
arc projection acquisition. A bidirectional encoder, which is used on the current gantry to 
track projection angle in the step mode, will be installed to track the projection angle on the 
arc. The tilt speed on the arc will be 7.57sce and the projection numbers on the arc will be 4 

25 U> 12. 
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. Since Ihc CIWCT system is based on an existing helical CT gantry and table, the 
system should have an existing computer-controlled table movement capability. With little 
modification, a circlc-plus-linc (CPL) scan can be achieved. Two bidirectional encoders are 
added: one is to track the longitudinal position of the x-ray source and the detector, and 
5 another to track the angular position of the source and detector. Then the system will be 
modified to synchronize x-ray exposure with 2 pulses/sec for line projection acquisition. 
Since only 25 cm of the trunk of the body in the z direction will be viewed per saan, the 
patient table is fed for 25 cm with the maximum feeding speed of 12.5 cm/sec, and then the 
volume scan time should be within 4-8 seconds, depending on the achievable feeding speed, 
10 required resolution and the actual size of the coverage in the z direction per scan. For 
detection of cancers such as lung cancer, circle-plus-lines and helical cone-beam scanning 
can also work. 

In addition to the method above to acquire circle and arc projections, alternatively, the 
circle-plus-arc geometry can be implemented in one of the following two ways. In the first 
15 and preferred of the three methods, the gantry 702 is tilted to a small angle (±15° to ±30°.) 
and then the x-ray tube 710 and the 2-D detector 71 1 are rotated while the gantry 702 is 
tilted. A half set of arc projections will be acquired only when the x-ray tube 710 and the 2-D 
detector 71 1 are at the rotation angles of 0° and 180°. When the tilted angle becomes zero, the 
circle projections will be acquired at the preset rotation angle positions. When the circle 
20 projection acquisition is completed, the gantry 702 will be tilted toward - 1 5° to -30°. Another 
halfsct of arc projections will be acquired only when the x-ray tube 710 and the 2-D detector 
71 1 are at the rotation angle of 0° and 1 80°. 

'Hie second alternative method is to mechanically modify a standard (T gantry such 
thai two short arc orbits arc added to Ihcr gantry, and the x-ray tube 710 mul the 2-1) detector 
25 71 1 can be moved on the arc to acquire the are projections and on the circle to acquire (he 
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circle projections. One nrc constitutes (he orbit of the x-ray lube 710 and the other are is the 
orbit of the 2-D detector 7 1 1 . The two arc orbits arc mounted 180° apurl from each other. The 
x-ray tube 710 and the 2-D detector 71 1 arc synchronously moved on the arc orbits to acquire 
arc projections. Then, the x-ray tube 710 and the 2-D detector 711 arc rotated on the gantry to 
acquire circle projections. 

Mounted on the gantry frame 702 opposite the x-ray source 710 is a 2-D detector 71 1 
which has a dynamic range equal to or greater than 1000:1 and an image lag of less than 
10%, for example a selenium thin film transistor (STFT) array or a silicon STFT array, in 
order to provide 2-D projections that correspond to an x-ray attenuation signal pattern. The x- 
ray source 710 and the 2-D detector 71 1 are mounted on the gantry frame 702 in such a 
manner that they both move synchronously. 

The cone-shaped beam of radiation 704 generated by the x-ray source 710 is projected 
through the body or objeot under test. The 2-D detector cone measures the radiation 
transmitted along the set of beam paths across the cone. 

Alternatively, a continuous series of two-dimensional detectors (not shown) can be 
fixedly mounted proximate to the gantry frame 702 and the x-ray source 710 is mounted to 
the gantry frame such that, upon rotation of the gantry frame, the cone-shaped radiation beam 
704 is projected through the body P under test and sequentially received by each of the series 
of detectors. 

A 2-D projection acquisition control and A/D conversion unit 726, under control of 
the scanning pulses sequentially obtained from the system control computer 724, which 
includes the clock 722, receives a sequence of outputs corresponding to different lines of the 
2-D detector 71 1 . Iiach line ofthe 2-D detector consists of many detection cells (at least 100). 
The output of each detector cell represents a line integral of attenuation values measurable 
liking one of the respective beam paths. The cone-shaped besim 704 subtends ti cone angle 
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sufficient to include the entire region of interest of the body. Thus, a complete scan of the 
object can be made by merely orbiting the gantry frame 702 supporting the x-ray source 710 
and the 2-D detector 71 1 around the body to acquire the 2-D projection signals at different 
angular positions. 

The analog-to-digital conversion unit 726 serves to digitize the projection signals and 
to save them in the 3-D image reconstruction array processor 728 and storage device 730. 
The method employed by the 3-D image reconstruction array processor 728 is the invented 
algorithm described herein. The 3-D image reconstruction array processor 728 serves to 
transform the digitized projection signals into x-ray attenuation data vectors. The x-ray 
attenuation data matrix corresponds to x-ray attenuation at spaced grid locations within the 
body trunk being examined. Each data element of the matrix represents an x-ray attenuation 
value and the location of the element corresponds to a respective 3-D grid location within the 
body. 

In accordance with the principles of the invention discussed previously, a display 
processor 732 obtains the data stored as 3-D x-ray attenuation signal patterns in the memory 
storage 730, processes the data as described above, and then the desired 3-D images are 
displayed on a 3-D display device 734. 

The 3-D image reconstruction array processor 732 may, for example, be comprised of 
an ULTRA SPARC- 10 model workstation, available from Sun Microsystems, Inc. of 
Mountain View, Calif. 94043. Another system is tho Mercury Computer Systems RACE 
Platform, which is a multiprocessor-based parallel computing system scalable up to a lew 
hundred processors. The reconstruction algorithm presented above is well suited to such 
parallel processing devices, since the various terms in the reconstruction can be calculated 
separately and summed. The use of a Storage Concept real-time storage system allows the 
acquisition ofup to 04 (315 ufdata continuously in real time. 
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The patient V is placed on n patient tabic 706 which is made to slide by a linear motor 

738 or some such device under control of the system control computer 724. Alternatively, 

the patient P can be placed on a fixed tabic, and a gantry frame holding the detector and the 

source can be moved over the patient P. 

An optional contrast solution injector 740, known in the art, can be used to inject a 

contrast solution for improved imaging. However, the invention can be used without such an 

injector. 

An example of the operation of the CBVCT tomography system 700 will now be 
explained with reference to Figs. 8-1 1. As shown in Figs. 8 and 9, the patient table 706 
bearing the patient P is moved into the gantry 702 so that the region of interest ROI lies 
between the source 710 and the detector 711. As shown in Fig. 10, to take the arc 
projections, the gantry 702 is tilted, and a cone beam 704 is emitted when the angular 
orientation of the source 710 is at 0° and 180° from a predetermined base location. As 
shown in Fig. 11, to take the circle projections, the gantry is righted, and the source 710 emits 
the cone beam 704 repeatedly as the gantry rotates. 

To decrease the total acquisition time, the sampling rate on the arcs can be reduced 
relative to the sampling rate on the circle. In addition, or as an alternative, the arc projections 
can be taken by using only a tilt-in of the gantry 702. A tilt-out of the gantry can be used to 
take additional arc projections to improve image quality. 

Developing and optimizing an x-ray scatter control and reduction technique is one big 
challenge for CBVCT because CBVCT is less immune to scatter than fan-beam CT. CBVCT 
image contrast is reduced by scatter without an effective control technique. Scatter can be 
countered with a hybrid technique that uses an air gap technique and an antimatter grid to 
control scatter and a practical software correction technique for detected scalier. One of the 
major differences between Ian beam slice (T and ClWCf is x-my beam collinuitioiu I Ising 
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very narrow slil collimntion in Ian beam CT reduces scutter-to-primary ratio (SPR) lo 0.2 or 
less. On live other hand, using a large oonc collimation in cone beam geometry with only an 
air gap technique results in an average SPR up lo 1. 

To overcome thai limitation, a software correction technique is used to correct for 
5 detected scatter and to reduce overall average SI'R lo 0.2 or less. Convolution filtering 
techniques and scatter detected by the FPD arc used to estimate scalier distribution and then 
subtract it from the total projection. A known convolution filtering technique taught in Love, 
LA, and Kruger, R.A., "Scalier estimation for a digital radiographic system using 
convolution filter," Med. Phys. 1987; 14(2): 178-1 85, was implemented for an image 
10 intensifier (Il)-based imaging system and produced an average percentage error of 6.6% for 
different anatomy and different clinical applications. That is equivalent to a reduction of SPR 
by a factor of up to 14. Even better scatter correction results can be achieved for an FPD- 
based system because there is no veiling glare component, compared lo an II-based system 
where that is a more dominant component. Based on previous studies and preliminary 
results, it is anticipated that the average SPR in each cone beam projection can be reduced to 
0.2. That is the equivalent SPR achievable in a fen beam slice CT, using a hybrid scatter 
correction technique (software correction plus air gap). That analysis and the preliminary 
results show that with the above-noted x-ray scatter reduction and correction techniques, the 
FPD-based CB VCTM system provides more than adequate low contrast resolution. 

The preferred embodiment combines an air gap technique wilh an antiscattcr grid and 
a software correction technique for residual scatter. A 10-15 cm air gap technique is an 
effective method to prevent large angle scalier radiation from reaching the dclcclor and to 
reduce average SPR to less than 2. It is contemplated that in the CHVCT system, the distance 
from the rotation center to the detector will he about 40 cm. With thai geometry, (he air gap 
is more (ban 15 cm lo achieve an average SPR less Hum 2. 
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One example of an efficient x-ray sealler rejection grid includes u focused, lunlalum, 
air-interspaced grid with a 10:1 grid ratio and 80 lines/inch. The grid strips are suspended 
between a pair of carbon fiber plates and aligned parallel to the axis of rotation. A scattcr-to- 
primary ratio (SPR) of approximately 1.0 can be achieved with 100 kVp and a moderate 
increase of the exposure level to keep the noise level unchanged. With a stationary grid there 
arc grid artifacts. To avoid such grid artifacts, the grid can be reciprocated with a computer- 
controllable speed to blur the grid strip artifacts. 

The residual scatter present within the projection images is removed based on a 
convolution-filtering method to estimate residual scatter distribution in each projection 
image. In the convolution filtering method, residual scatter is modeled as a low pass, 
spatially filtered version of the total projection (scatter plus primary). After estimating 
residual scatter in each projection, the residual scatter radiation is then subtracted to obtain 
primary distribution for reconstruction. That technique effectively reduces SPR from 1.0 to 
0.2 or less. 

The conventional convolution filtering method requires two x-ray projections at each 
projection angle to accurately estimate residual scatter: one with a beam stop array for 
calculating two scaling factors and another without the beam stop array. That is not practical 
and would significantly increase patient dose in CBVCT. To overcome those difficulties, the 
preferred embodiment uses scout images for estimating scatter distribution in "real time" for 
each patient. Before starting to scan, one scout projection image is acquired, as in a standard 
fan beam CT. Traditionally, the scout images are used for positioning, and surveying body 
size to adjust the x-ray exposure levels in real time and reduce patient dose Before acquiring 
scout images, as Shown in Figs. 12A and I2B, a square matrix 1204 of small lead ball 
bearings 1200 is placed between the x-ray collimator 1202 and the region of interest KOI. 
Both primary and sampled scalier distributions are estimated from the scout images with the 
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lead beam stop array. The estimated primary images are used lor a scouting purpose. The 
scaling factors lor estimating scatter distribution and the convolution kernels at sampled 
angle positions can be determined. Then the scatter distributions are estimated using the 
convolution kernel at corresponding angle positions and subtracted from the detected 
projections. To reduce radiation dose to the patient and computation load, only a minimum 
number of required scout images are acquired. Only a few scout images are needed because 
the accuracy of the method is not highly dependent on the exact shape of the convolution 
kernel, so long as its dimensions are large enough. The exponential kernel is used for the 
estimation of residual scatter because a 2D exponential kernel is an optimum formation. 

Another technique which can be used in the present invention to improve imaging is 
the ultra-high-resolution volume-of-interest (VOI) reconstruction mode. That technique can 
be used to focus on a suspicious lesion. 

It is known in the art for fiat panel detectors to have zoom modes. One source of such 
flat panel detector is Varian Imaging Products of Mountain View, California, U.S.A. The 
Varian PaxScan 2520 flat panel detector has the following characteristics: size = 19.5 x 24.4 
cm, frame rate = 15-120 frames per second, image lag < 10°/o, pixel pitch ~ 127 fim, A/D « 
16 bits, exposure range = 1-3000 uR, DQE = 65%, dynamic range = 2000-30,000:1. Even 
larger flat-panel detectors are known in the art, e.g., 50 cm x 50 cm. 

The zoom mode of a flat panel detector such as a Varian flat panel detector is used to 
acquire projection data for ultra-high VOI reconstruction. In the zoom mode, the detector can 
acquire a random block of 768 x 960 pixels at 30 framcs/scc. with the full 4 Ip/mm resolution 
of the sensor. The pixel size of the detector, as noted above, is 127 [tin. A dual-focus spot x- 
ray tube is used, having focus spots of 0.3 and 0.6 mm. Ullra-high-resolulion VOI can use a 
0,3mm focus spot, so that the focus spot size will not be a limiting factor of the spatial 
resolution for the VOI mode. Therefore, the I'OV (Held of view) oflhe zoom mode is 9.75 x 
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12.2 cm. To reduce unnecessary radiation (o (lie patient, a collimator limits the radiation to 
within the ROI (region of interest) in the VOl acquisition. A narrow strip of collimation (-2 
cm wide) is needed. If the ROI is larger than 12.2 cm in diameter, the projection data 
acquired in ultra-high VOI mode arc truncated in the lateral direction. There are some streak 
5 artifacts if the reconstruction is obtained from the truncated data without preprocessing the 
data. The conventional method to deal with truncated projection data is to tail the projection 
data with a cosine wave before filtering (Z. H. Cho, E. X. Wu, S. K. Hilal: "Weighted 
backprojection approach to cone-beam 3D projection reconstruction for truncated spherical 
detection geometry," IEEE Trans Med Imaging 13(1), 1 10-122, March, 1994). Fortunately, 

10 in the present case, the complete information in the region out of VOI is already available 
from the previous lower resolution scan. That information can be used to tail the truncated 
projection data and then complete the VOI reconstruction. Computer simulation indicates 
that such an algorithm eliminates the reconstruction artifacts introduced by truncated data 
within VOI. Such a technique is anticipated to be better than the conventional method. It is 

15 further anticipated that the ultra-high-resolution VOI reconstruction technique can provide up 
to 5.0 Ip/tnm resolution with a justifiable increase of the x-ray dose. The above-disclosed 
VOI technique can be used to detect cancers, such as breast and lung cancer. 

A FPD-based CBVCT system will provide better contrast and spatial resolution and 
better geometric accuracy than an II-ba$ed CBVCT system. Recently, a new technology of 

20 large area flat panel solid state detector array has been developed by several groups. A high 
resolution, high frame rate, amorphous silicon (a-Si:Il) FPD using a phosphor screen and a 
photodiodc array to convert incident x-rays to a charge image has been developed. Also, a 
selenium VPD has been developed by other groups using a uniform layer ofan x-ray sensitive 
pholoeoiuluctor, selenium, for a direct conversion of x-rays to an electronic image, in 

25 addition, a real time WD for fluoroscopic images has been developed. In spite of their 
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differences, these image sensors have some common potential advantages over other 
detectors: compactness, high DQB. absence ol" geometric distortion and veiling glare with the 
benefits or high resolution, high frame rate, high dynamic range, small imngc.lag (< |%) an d 
excellent linearity (-1%). The F1>D has almost the same DQB as tin II within the diagnostic 
radiation range. These advantages of the new FPD over an ll-CCD detector make it a good 
candidate for the detector used in CBVCT. Therefore, a FPD-based CBVCT system will 
make CBVCT a superior technique. In the past two years, the development of the TFT 
detector has been exciting and progressed from the research phase to the production phase. 
Six companies have started to manufacture this type of detector. 

The FPD-based CBVCT angiography system has better spatial resolution and low 
contrast resolution than Il-bascd systems. The FPD-based system has better spatial resolution 
than a helical CT and near equal low contrast detectability in comparison to a helical CT. 

Two alternatives to the circle-plus-multiple-arcs orbit can be used to obtain data 
sufficient for exact reconstruction. 

Fig. 4C shows an orbit for taking a scan over a circle plus multiple (in this case, two) 
lines. That orbit allows at least 25 cm coverage in the Z direction. Before starting the scan, 
the patient on the table is positioned to -12.5 cm from the center of the scanner to prepare for 
circle-plus-line (CPL) scan. While the table is moving toward the center of the scanner, the 
x-ray tube-detector rotates, taking projections only at 0° (on the upper line) and 180° (on the 
lower line) or the rotation angle positions to obtain two line projections per rotation. When 
the table is at the center or the scanner, the toblc stops moving, and the x-ray tube and 
detector rotate to acquire multiple circle projections. Alkr acquisition or (he circle 
projections, the table moves toward the position of + 12.5 cm from the center of the scanner, 
while the x-ray tube and detector are rotating, taking more line projections as ahove. The 
reconstruction algorithms are those taught by I In. 
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To reduce CTL scan time, Ihc "quasi" spiral scan mode of lite gantry is used because 
during the scan, Ihc x-ray tube and detector continue rotating on the gantry while the table is 
moving and the table slops at the center of the scanner to acquire circle projections. This 
scan mode will eliminate x-ray lube-detector rotation-stops during scan and reduce (he 
transition lime between line acquisition and circle acquisition. II can also be noted that we 
actually acquire line projections along two lines: the upper and lower lines, as shown in 
Figure 4C which shows the implementation of CPL orbit using "quasi" spiral scan on a spiral 
gantry with 8 line projections at the positions numbered 1 through 8. This is because this 
can reduce the sampling rate on a single line and increase line-scanning speed when using 
"quasi" spiral mode. 

A spiral scan will now be explained with reference lo Fig. 4D. While the table is 
moving through the gantry, the x-ray tube and detector rotate, talcing projections at multiple 
angle positions to obtain multiple spiral cone beam projections per rotation. Depending on the 
achievable frame rate of the detector, the spiral scanning time will be 2-8 seconds and the 
total spiral rotation angle will be from 1 80° plus cone angle to 720° to cover 6.5 cm to 25 cm 
in the z-direction. The reconstruction algorithms are those taught by Wang, G.E., Lin, T.H., 
Chen, P.C., and Shinozaki, D.M., "A General Cone-Beam Reconstruction Algorithm/* IEEE 
Transactions on Medical Imaging, Vol. 12(3):486-496 (1993), and in Wang, O.E., Lin, T.H., 
and Chen, P.C., "Half-Scan Cone-Beam X-Ray Microtomography Formula," Scanning, Vol. 
16, pp. 216-220(1994). 

Based on currently available FPDs, a spiral scan should have a scanning speed of 2-8 
seconds, a z-covcrage of 65-250 mm and a slice thickness in the /. direction of 0. 1 7-0.67 mm. 
That scan offers the advantages of more uniform sampling and ease of implementation. Hie 
circle-plus-arcs and circle-plus-lines scans should both have a scanning speed of 4-8 seconds, 
a /-coverage of 130-250 mm and the same slice thickness as thai for the spiral scan. The 

38 
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circic-plus-liiics ori^it is currently the best to nddress Iruncution. The circle-plus-urcs orbit 
has the advantage that it requires no patient transition during the scan; that is, the patient 
remains still. 

The scanning limes listed above arc estimated based on the frame rate of currently 
available FPDs (60 -120 framcs/s) and the spiral gantry speed for a used spiral CT gantry (I 
s/rcvolution). Existing FPDs arc specially designed for radiographic or fluoroscopic imaging 
not for fast tomographic imaging. Once a large size FPD-based CBVCT becomes feasible, a 
large size FPD specially designed for fast tomographic imaging with high speed and low 
image lag will be developed. If the frame rate of the detector is increased to 960 frames/sec. 
for a 25 cm x 50 cm FPD, the spiral scan time for 33 cm coverage will be I - 2 seconds 
depending on gantry speed and table moving speed. 

To deal with the projection truncation problem, the following three measures can be 
taken. First, when determining which orbit and . related algorithm will be used for 
reconstruction, the one with the smallest contaminated depth should be used. Second, when 
measuring the total dose to a patient, the dose received in the contaminated region due to 
projection truncation should be included. Third, when acquiring projection data using the 
CPL or CPA, the detector in the longitudinal direction should be slightly larger than the ROI. 

If it is necessary for the ultra-fast readout of FPD, a subtraction algorithm can be used 
to reduce the effect of image lag. In such an algorithm, the previous N weighted projections 
will be subtracted from the current projection. The weighting factor for each previous 
projection will be determined by the lag measured vs. frame numbers subsequent to the frame 
in which it was generated. Then the final image will be reconstructed from image lag- 
corrected projections. 

While a preferred embodiment of the present invention has been set forth above, those 
skilled in (he art will recognize thai other embodiments can l>e realized wilhin the scope of 



V) 



WO 02/30282 PCT/US01/32011 

the invention. Tor example, numerical values and the names of speeifie products are 
illustrative rather than limiting. Also, to achieve the circle-plus-two-arcs or circlc-plus- 
multiplc-arcs orbit, any suitable equipment and mode of operating it can be used. 
Furthermore, the particular algorithms presented herein are illustrative rather than limiting. 
Moreover, while the utility of the invention has been presented with particular attention to 
lung cancer, it can be used for other malignancies. Therefore, the invention should be 
construed as limited only by the appended claims. 
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AIM'KNDIX 

A. Derivation of the Reconstruction Algorithm for Arc CB Projections 

( x 4*y«'**) is lhc fixed coordinate system, and (x,o>,,z,) is the local coordinate 
5 system rotating rigidly in phase with the virtual detector plane. The arc orbits arc within the 
plane determined by (x 0 ,y lt ). In the local coordinate system (x t% y L ,z f ), Grangcafs 
formula can be written as: 

I O^i I 1 c) 2 * * l I CY) I A 

- 1 « { y^P»0\Z$0'smQ + Z C osQ-l)dYdZ = ^Rf(OSn,n) (A-I) 

where Rf(OS h,n) is the Radon transform over the shadowed plane with norm vector n 
10 which passes through the source focal point S, ]',and Z t are the integral limits along the J' and 
Z axes respectively, and P p ()\Z) is the cone beam projection at angle p along the arc 



orbits. 



In the local coordinate system (x^y,^) , 



OS = (A 0, 0) 



(A-2) 
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n as (sinG cos(<p - p),sin0 sin(cp - p),cos0 ) 



(A-3) 



p = OS • n = £>sin0 cos(<p - p ) 



(A-4) 



|a9xw| 2 =D a -D I sin 2 Gcos i (<p-p) = /) 2 -p 2 



(A-5) 



OS -n = r*n 



(A-6) 



For derivation convenience, by letting 
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(A-I) is converted into 
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(A-8) 



Taking the I s ' derivative on both sides of(A-8) and using (A- 1 5) (see below) gives 



D 1 



iy--p> 



2p 0 



/ ( (p,/,0) = iL./{/(rn.r7) (A-9) 



In the local coordinate system {x L ,y L ,z L ), the Radon plane SD t Dj can be 
5 represented by the equation: 

lx L + DcosBy, + £>sin ©r, - D/ = 0 (A-10) 
On the other hand, in the fixed original coordinate (x (t ,y 0 ,z 0 ), another equation to describe 
the same Radon plane can be written as 

x H sin9 costp + y () sinG cos(p + z 0 cos0 - p = 0 
10 Through the relation between the local and the fixed coordinate system 
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and solving (A-10) and (A-ll) simultaneously, we have: 
- D cos 0 sin p -f / cos (i 
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sin9 coscp = 

sin0 sincp = 
Dl 



Dcos0cosp 4-/sin ft 



cost) 



D sin0 



Consequently, I'rom (A- 1 3) and (A-M), we have 



(A-13) 
(A-M) 
(A-15) 
(A-16) 
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0 cos 6 con It +/sinB 
tunro = ! !— 

-0cos0Kin(H-/cos|l v ' 

Considering the variable change from (0,<j>) (o.fl) , we gel 
. . _ . /Jcos0</0 m 

(g'w'l^ (A " 18) 

Furtlicr, by introducing (A- 15) into (A-9), wc get 

^^[^ + -^-^jA(P./.e) •„,*) (A . 19) 

According to the 3D Radon inverse transform (27), eventually, we have 

l v'f 5 2 

= "4^ J J ' r, ) sin0de</( P (A-20) 

where w{pi/ s 0) is the window function for the sub-domain missed by the circular orbit in 
the Radon domain, / = {1,2} corresponds to arc orbit 1 and arc orbit 2 respectively, and p, 
determines the arc sub-orbit spanning range. 
B. Derivation of the Window Function 

By referring to Fig. 3, the section in Radon domain that can not be covered by the 
radon transform of the circular CB projections is 

Q * a |",p]|/Xl-sin 2 0siir<p)^ <|p| 
Actually, we have (see Appendix A) 

/;/ 

p= 7/7-m-/-- < ,J - 2 > 
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cosO 



tancp = 



DcosQcosP +/sinP 



- DcosGsin p +/cosp 
Then, from (B-3) We gel 



sin z 0 = l-cos 2 0 = r — ; — 



5 and from (B-4) we get 

• 1 te ? q> _ QPcos©cos p +/sin p ) 2 
Sm 9 ~l + lg\ ~ D 3 cos 2 0+/ 2 

Incorporating (B-2), (B-5) and (B-6) in (B-l) we get 



=>Z> 



1- 



(Z) cos© cos p +/sinp) 2 
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Z) 7 4/ 2 

|/|> [d 2 +/ j -(DcosQcosp +/sinp) 3 )^ 

/ 2 >D J +/ 2 -(Z)cos0cosP +/sinp) 2 

D 2 <(£>cos©cosP +/sinp) 2 
Finally, the window function can be written as 



«<P,/,0)= 



1 /sin P >£)•(! -cos0cosP) 
1 /sinp <-£» 0 + cos©cosP) 
0 elsewhere 
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I claim: 

I. A method of imaging an object to form a reconstructed image, the method 
comprising: 

(a) scanning the object using a source of radiation and a detector or the radiation, the 
source and the detector being moved around the object lo define an orbit 
comprising (i) a circle orbit for providing a first set of data signals and (ii) a 
plurality of arc orbits for providing a second set of data signals; 

(b) performing a first reconstruction from the first set of data signals lo generate a 
first reconstruction result; 

(c) performing a second reconstruction from the second set of data signals to generate 
a second reconstruction result; and 

(d) summing the first reconstruction result and the second reconstruction result to 
obtain the reconstructed image as a sum of the first reconstruction result and the 
second reconstruction result 

2. The method of claim I, wherein the source of radiation is a source of cone- 
beam radiation. 

3. The method of claim 2, wherein the first reconstruction result is a sum of a 
Feldkamp reconstruction and a complementary term used to correct the Feldkamp 
reconstruction for use of the cone-beam radiation on the circle. 

4. The method of claim 2, wherein the second reconstruction result is generated by 
summing individual arc reconstruction results for the arc orbits. 

5. The method of claim 4, wherein two arc orbits arc used, 

6. The method of claim 5, wherein, for each or the arc orbits, the individual arc 
reconstruction result is culcululcd using a window function to account for a 
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portion of a Radon domain which cannot be covered by a Radon transform of the 
circle orbit, 

7. The method of claim 6, wherein, for any point r in the object: 

(*/, J'£. zi) is a detector coordinate system which rotates rigidly with the detector; 
P is an angle along the arc orbits at which the source is located; 
D is a radius of the arc orbits; 

r and a location of the source define a plane which intersects a plane of the 
detector at a line of intersection, the line of intersection having a point C of closest 
approach to an origin O of the detector coordinate system such that a line segment 
connecting O and C has a length / and forms an angle 0 with a y L axis of the 
detector coordinate system; and 
the window function is given by 



w(p./ f e)« 



I /sinP > £)-(l-cos0cos0) 
1 /sinp <-Z)(l + cos©cosP). 
0 elsewhere 



8. The method of claim 7, wherein: 

{Y t Z) is a projection of r onto the detector coordinate system; 

Y f and Z, are integration limits determined by dimensions of the detector; 

P, is an integration limit determined by dimensions of the arc orbits; 

/.•(P,A8)« J // > n . l (K,Z)5(rsjn© + Zcos0-/>/raZ; 

,.,,K,e) = ^,,( p .,e,co s9 [|.| t i^|L],, ([ ,,e ki „, d 

the individual arc axons! ruction results arc given by 
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{ )> The method ofclaim X, wherein: 

Mi) is an integration limit determined by a spatial sampling frequency of the 
detector; 

Ly m is an integration limit along a Z direction; 

/\t»( J'Z) is a conc-bcam projection at a point (Y,Z) on the detector; 

M z )- Jh|^exp(/a)Z); 





tiie Feldkamp reconstruction is given by 




-/>,;(y(r),Z(F)); 



A y»(j')= }yo>cxp(/<»r)rf<o ; 



^•(i')=^(i'.*>ft; 



-/.. 



the complementary term is given by 



,(^) = -^r <f cAIi— 




the first reconstruction result is given by 
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10. The method of claim 9, wherein step (b) is performed through a digital signal 
processing operation which assumes: 



I 

4 • " = 0 

'-(-»)" I W ' ,h «*0 

87T 1 » 3 



and 



° I •„ " = ° 



11. The method of claim 1, wherein the detector is a flat panel detector. 

12. The method of claim 1, wherein at least one of steps (b) and (c) comprises 
scattering correction. 

13. The method of claim 1, wherein the first reconstruction result and the second 
10 reconstruction result are both in a Filtered backprojection format. 

14. The method of claim 1, wherein steps (b)-(d) are performed through parallel cone 
beam reconstruction. 

15. The method of claim 1, wherein the object is longitudinally unbounded, and 
wherein step (d) comprises providing the reconstructed image as an exact 

15 reconstruction of the longitudinally unbounded object. 

16. The method of claim 1, wherein step (d) comprises forming a 3D matrix of 
attenuation coefficient distribution of the object. 

17. The mclhod of claim I, wherein the object is a patient or a region of interest in a 
patient. 

20 1 8. The method of claim I, wherein the reconstructed image is formed for 

nondestructive testing of the object. 
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19. The method ofclnim I, whercin step (a) comprises performing a qtmsi-spiml scan 
olthe object. 

20. 'l-hc method of claim 19, wherein the quasi-spiral scan comprises a tilt in p | us 
circle scan. 

21. The method olclaim 20, wherein the quasi-spiral scan further comprises a tilt out 
scan. 

22. The method of claim 1, wherein the first set of data signals has a higher sampling 
rate than the second set of data signals. 

23. The method of claim 1, further comprising: 

(e) locating a volume of interest in the object; 

(f) scanning the volume of interest; and 

(g) performing steps (b)-(d) again for the volume of interest. 

24. The method of claim 23, wherein the detector has a first resolution and a 
- second resolution higher than the first resolution, and wherein the first resolution is 

used in step (a) and the second resolution is used in step (e). 

25. A device for imaging an object to form a reconstructed image, the device 
comprising: 

a source of radiation; 

a detector of the radiation; 

a gantry for supporting the source and the detector and for moving the source and 
the detector around (he object to define an orbit comprising (i) a circle orbit for 
providing a first set of data signals and (ii) a plurality ofarc orbits for providing a 
second set ofdatu signals; and 

a computing device, receiving the first and second sols of data signals, for (i) 
performing :i first reconstruction from the first set of data signals to generate a first 
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reconstruction result, (ii) performing a second reconstruction from the second set of 
data signals to generate a second reconstruction result and (iii) summing the first 
reconstruction result and llic second reconstruction result to obtain the reconstructed 
image as a sum of the first reconstruction result and the second reconstruction result. 

26. The device of claim 25, wherein the source of radiation is a source of conc-bcam 
radiation. 

27. The device of claim 26, wherein the first reconstruction result is a sum of a 
Feldkamp reconstruction and a complementary term used to correct the Feldkamp 
reconstruction for use of the cone-beam radiation on the circle. 

28. The device of claim 26, wherein the second reconstruction result is generated by 
summing individual arc reconstruction results for the arc orbits. 

29. The device of claim 28, wherein two arc orbits are used. 

30. The device of claim 29, wherein, for each of the arc orbits, the individual arc 
reconstruction result is calculated using a window function to account for a portion of a 
Radon domain which cannot be covered by a Radon transform of the circle orbit. 

3 1 . The device of claim 30, wherein, for any point r in the object: 

Xxi* yh zi.) is a detector coordinate system which rotates rigidly with the detector, 
p is an angle along the arc orbits at which the source is located; 
D is a radius of the arc orbits; 

F and a location of the source define a plane which intersects a plane of the 
detector at a line of intersection, the lino of intersection having a point C of closest 
approach to an origin O of the detector coordinate system such that a line segment 
connecting O and C has a length / and forms an angle 0 with a y tt axis of the 
detector coordinate system; and 
the window function is given by 
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\ /sin fl > O (l-cos0cos|l) 
I /sin <«O (l+cos0cosH). 
0 elsewhere 



32. The device of claim 31 , wherein: 

•(}', Z) is a projection of r onto the detector coordinate system; 

Y, and Z f arc integration limits determined by dimensions of the detector; 

(it is an integration limit determined by dimensions of the arc orbits; 

/,(P,/,0)= f J/ 5 ^ l (> / >2:)5(rsin0 + Zcos0-./>/}yZ; 

,J P xe).^l,Kp, ( ,9) MS e[|.| + ^L|.], (p ,,, e);ai , d 



the individual arc reconstruction results are given by 
I 



33. The device of olaim 32, wherein: 

0) D is an integration limit determined by a spatial sampling frequency of the 
detector, 

Lz is an integration limit along a Z direction; 

PutY.Z) is a cone-beam projection at a point (Y,Z) on the detector; 
to, 

Jl»|*)cxp(/<DZ); 

'C(»'.z)= Jifc/i H (^-=)/;(r,r); 

(he Kcklk.-imp reconstruction is given by 
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a, 

>'mO')= Jycocxp(/<D I'Vo); 

-/., 

5 the complementary term is given by 

Aft— TJT f * V n y^ W): and 
471 [o.i„l (£>-r*J 

the first reconstruction result is given by 

34. The device of claim 33, wherein the computing device generates f t (r) through i 
1 0 digital signal processing operation which assumes: 



fy,|(")= 



with 



/i = 0 



4 

'-(-!)" 1 **0 
8n 2 /» 2 



and 



H)" 1 w,lh **0 

/J 



35. The device of claim 25, wherein the detector comprises a Flat-panel detector. 
15 36. The device of claim 35, wherein the flat-panel detector has a frame rate of at least 30 
frames per second, 

37. The device of claim 35, wherein the Hat-panel detector has an image lag of less than 2%. 
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38. The device of claim 35. wherein (lie Hal-panel detector lius u dynamic range of ut least 
4000:1. 

39. The device of claim 35, wherein the flat-panel detector has a size of at least 19,5 cm x. 
24.4. cm. 

5 40. The device of claim 25, further comprising a device for correcting for scatter in the 
reconstructed image. 

41. The device of claim 40, wherein the device for correcting scatter comprises an antiscatter 
grid disposed between the source and the object. 

42. The device of claim 41, wherein the device for correcting scatter further comprises a 
10 collimator for collimating the radiation. 
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Figure 1 
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Figure 8 
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